Non-LTE line formation for Fe I and Fe II 
in the atmospheres of A-F type stars 



L. Mashonkina 

Institute of Astronomy, Russian Academy of Sciences, RU-119017 Moscow, Russia 

o 

(N 

I Abstract. Non-local thermodynamical equilibrium (non-LTE) line formation 

04 ■ for neutral and singly-ionized iron is considered through a range of spectral types 

when the effective temperature varies from 6500 K up to 8500 K, the gravity 
from log g = 4 down to log g = 3, and the metallicity is solar one. The non-LTE 
calculations were performed with a comprehensive model atom for iron which was 
treated in our earlier paper (Mashonkina et al. 2011). The departures from LTE 
lead to systematically depleted total absorption in the Fe I lines and positive 
, abundance corrections, in qualitative agreement with the Rentzsch-Holm (1996) 

%^ \ study. However, our predicted magnitude of the non-LTE effects is significantly 

c/3 ■ smaller compared to the previous results due to the use of rather complete model 

^ atom of Fe L The non-LTE abundance corrections do not exceed 0.1 dex for the 

dwarf models and 0.20 dex for the giant ones. Non-LTE leads to strengthening 
^ ; the Fe II lines, however, the effect is small, such that the abundance correction 

■ is at the level of —0.01 to —0.03 dex over the whole range of stellar parameters 

^ ■ being considered. No firm conclusion can be drawn with respect to whether or not 

t:^- ! the Fe I/Fe II ionization equilibrium is fulfilled in the atmosphere of the Sun and 

I Procyon due to uncertainty in available gf— values for visible lines of Fe I and 

O ■ Fe II. 
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1. Introduction 

Iron is the most popular chemical element in studies of both non-magnetic and magnetic A and later 
type stars thanks to quite numerous lines of the neutral and singly-ionized species in the visible 
spectrum. Iron lines are used to determine stellar parameters, such as effective temperature, Tefj, 
surface gravity, log g, microturbulence velocity, Vmic, and, in the magnetic stars, also to investigate 
magnetic fields. Therefore, the calculations of the iron lines have to be based on realistic line 
formation scenarios. 

In the atmospheres with T^s > 4500 K, Fe I is a minority species. Number density of the minority 
species can easily deviate from the thermodynamical equilibrium (TE) population due to small 
deviation of the mean intensity of ionizing radiation from the Planck function. Since the beginning 
of 1970th a number of studies attacked a problem of the non-local thermodynamic equilibrium (non- 
LTE) line formation for Fe I in the Sun and late type stars (Tanaka, 1971; Athay & Lites, 1972; 
Boyarchuk et al., 1985; Takeda, 1991; Gratton et al., 1999; Thevenin &: Idiart, 1999; Gehren et al., 
2001; Shchukina & Trujillo Bueno, 2001; Collet et al., 2005) Here we list only the studies, where 
original model atom was produced. The non-LTE effects in the hotter atmospheres were investigated 
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by only Gigas (1986) and Rentzsch-Holm (1996). All these research came to a common conclusion 
that Fe I is subject to overionization in stellar atmospheres due to enhanced photoionization of the 
levels with the ionization edges in the ultra-violet (UV) and the Fe I lines are weakened compared 
to their LTE strengths. However, a consensus on the expected magnitude of the non-LTE effects 
was not reached. For example, the surface gravity correction for cool metal-poor atmospheres was 
calculated, in some studies, at the level of 0.5 dex, while no significant departures from LTE were 
found by others. The discrepancies between different papers might be due to an incompleteness of 
the applied model atoms. 

Mashonkina et al. (2011) constructed a comprehensive model atom for iron using more than 3000 
measured and predicted energy levels and showed that the departures from LTE for Fe I decreased 
significantly in the atmospheres of solar-type stars, in particular, in the metal-poor atmospheres 
compared to those calculated in the previous non-LTE studies. In this study, we revise the non- 
LTE effects for the stellar parameter range characteristic of A-type stars. The paper is organized 
as follows. An introduction to the non-LTE approach is given in the next section. Sect. [3] describes 
briefly the method of the non-LTE calculations for Fe I-Fe II. The results for the two reference 
stars, the Sun and Procyon, are presented in Sect. HI Finally, we report on the departures from LTE 
depending on stellar parameters. 

2. What is meant by a non-LTE approach? 

We assume that the particle velocity distribution is Maxwellian with a single kinetic temperature 
for various kind particles, i.e., for electrons, atoms, and ions: Tg = = Tj. We consider that the 
occupation number of any atomic level is determined from the balance between radiative and colli- 
sional population and de-population processes, such as photoexcitation, photoionization and their 
inverse processes, inelastic collisions with electrons, atoms, molecules, dielectronic recombination, 
charge exchange, i.e., from a statistical equilibrium (SE): 

niY,{R^■j + Cij) = Y.n,{R,, + Cji) i = l,...,NL. (1) 

Collision rate Cij is defined by the local kinetic temperature Tg and collider particle number density. 
Collision processes tend to establish TE. Radiative rate Rij depends on the radiation field, which 
is, in general, highly non-local in character. Radiative processes tend to destroy TE. 

In the nature, each chemical species possesses a huge number of bound energy levels approaching 
to the infinity. In practice, the real atomic term structure is represented by the model atom with 
finite number of levels NL. To investigate the SE of an atom, a large amount of various atomic 
data is required, such as energy levels, photoionization cross-sections for every level in the model 
atom, oscillator strengths (fij) for the entire set of allowed transitions, and collision excitation and 
ionization data for the entire set of transitions. For the last two decades, the situation with atomic 
data has improved significantly. A fairly extensive set of accurate atomic data on photoionization 
cross-sections and oscillator strengths was calculated in the Opacity Project (Seaton et al., 1994) 
and has become available via the TOPBASE database. The IRON project is in progress and provides 
the data for the Fe group elements. Quantum-mechanical ab initio calculations exist for electron 
impact excitation in a number of selected atoms and ions. However, for heavy elements beyond the 
Fe group, much data are missing, in particular, for the high-excitation states. Rough theoretical 
approximations are still used to evaluate inelastic collision cross-sections. The calculated SE of an 
atom depends on the completeness of the adopted model atom and the accuracy of the used atomic 
data. 

The excitation and ionization states of the matter is calculated from the solution of combined 
SE ([T]) and radiation transfer equations ([2]): 
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M^^^^ = -X.(^)/.(^,/x) + Mz) (2) 

Here, Xi^ and rji, are the absorption and emission coefHcients, respectively. The radiation transfer 
equations have to be solved at frequences of all the radiative transitions in the atom which makes the 
non-LTE calculations to be bulky. In non-LTE, atomic level population ni{d) at any depth point d de- 
pends on the physical conditions throughout the atmosphere: ni{d) = f{ni, . . . ,nNL, Ji, ■ ■ ■ , Jnf)- 
Here, J,y is the mean intensity. The models in which the Saha-Boltzmann equations are replaced by 
the physically more accurate SE equations are called the non-LTE models. 

For comparison, in LTE, level populations are calculated from the Saha and Boltzmann equations 
and depend only on local temperature and electron number density. This implies simple numerical 
calculations and the need in atomic data only for the spectral line under investigation. The LTE 
assumption is valid if every transition in the atom is in detailed balance. This is not fulfilled in 
spectral line formation layers, where the mean- free path of photons is large and the intensity of 
radiation is far from being in TE. To evaluate the effect of the departures from LTE on line strengths 
and stellar parameters derived from these lines, one needs to solve the non-LTE problem for a given 
chemical species. 



3. The method of non-LTE calculations for iron 

The non-LTE calculations were performed using a revised version of the DETAIL program (Butler 
& Giddings, 1985) based on the accelerated lambda iteration (ALI) scheme. 

Model atom of Fe I-Fe II. We applied the model atom treated by Mashonkina et al. (2011). It was 
constructed using the 958 energy levels of Fe I known the experimental analysis of Nave et al. (1994) 
and also the levels predicted from calculations of the Fe I atomic structure by Kurucz (2009), in 
total 2970 levels. The system of measured levels is nearly complete below £^exc — 5.6 eV. However, 
laboratory experiments do not see most of the high-excitation levels with i?cxc > 7.1 eV, which 
should contribute a lot to provide close collisional coupling of Fe I levels to the Fe II ground state. 
Neglecting the multiplet fine structure gives 233 levels of Fe I in the model atom. The predicted 
high-excitation levels were used to make up six super-levels. This is new compared to all previous 
non-LTE studies for Fe I. In total, 11958 allowed transitions occur in this model atom of Fe I. For 
Fe II, we use the 89 lowest terms with E'exc up to 10 eV. The ground state of Fe HI completes the 
system of levels in the model atom. 

Photoionization is the most important process deciding whether the Fe I atom tends to depart 
from LTE in the atmosphere of a cool star. Our non-LTE calculations rely on the photoionization 
cross-sections of the IRON project (Bautista ct al., 1997). The hydrogcnic approximation was used 
for the Fe II states. Radiative bound-bound (b-b) transition rates were computed using (^/-values 
from the Nave et al. (1994) compilation and Kurucz (2009) calculations. All levels in our model 
atom are coupled via collisional excitation and ionization by electrons. The collisional rates were 
calculated using the theoretical approximation of van Regemorter (1962) and Seaton (1962). 

Model atmospheres. For the Sun and Procyon, the calculations were performed with plane-parallel 
and blanketed LTE model atmospheres computed with the MAFAGS-OS code (Grupp et al., 2009), 
which is based on up-to-date continuous opacities and includes the effects of line-blanketing by means 

of opacity sampling. Small grid of model atmospheres with T^s = 6500 — 8500 K and log g = 3 and 
4 was taken from the Kurucz (2008) website to evaluate the departures from LTE depending on 
stellar parameters. 
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Figure 1: Departure coefficients, b, for the levels of Fe I and Fe II as a function of logrsooo in the 
model atmosphere 8000/4/0 from the calculations with our final model atom (left panel) and the 
reduced model atom which ignores the predicted levels of Fe I (right panel) . 



Statistical equilibrium of iron. The departure coefficients, bi = n^^'^^/nY^^, for the levels of Fe I 
and Fe II in the representative model atmosphere T^g/log g/[M/H] = 8000/4/0 are presented in 
Fig. [1] as a function of logrsooo- Here, ^^^lte ^^^^ nY^^ are the statistical equilibrium and thermal 
(Saha-Boltzmann) number densities, respectively. The difference between two panels of Fig. [T] is in 
using two different atomic models, i.e., the final model atom from Mashonkina et al. (2011) in the 
left panel and the reduced model atom which ignores the predicted levels of Fe I in the right panel. 
Our calculations support qualitatively the previous results of Rentzsch-Holm (1996) for common 
stellar parameters. 

• All the levels of Fe I are underpopulated in the atmospheric layers above log T5000 = —0.3 due 
to the overionization caused by superthermal radiation of a non-local origin below the thresholds of 
the low excitation levels of Fe I. 

• Fe II dominates the element number density over all atmospheric depths. Thus, no process 
seems to affect the Fe II ground-state and low-excitation-level populations significantly, and they 
keep their thermodynamic equilibrium values. 

In this study, progress was made in establishing close collisional coupling the Fe I levels near the 
continuum to the ground state of Fe II, due to including in the model atom a bulk of the predicted 
high-excitation levels of Fe I. As a result, the populations of the Fe I levels obtained with the final 
model atom (Fig.[Tl left panel) are closer to the corresponding TE populations in the line- formation 
layers (logrsooo = up to —3) than that for the reduced model atom (Fig.[T] right panel). 

4. Fe I/Fe II ionization equilibrium in the Sun and Procyon 

The non-LTE method was applied to analyze the Fe I and Fe II lines in the two stars with well- 
determined stellar parameters: the Sun and Procyon. The solar flux observations are taken from 
the Kitt Peak Solar Atlas (Kurucz et al., 1984). The spectroscopic observations for Procyon were 
carried out by Korn (2003) with the flbre-fed echelle spectrograph FOCES at the 2.2m telescope of 
the Calar Alto Observatory, with a spectral resolving power of ii ~ 60 000 and a signal-to-noise of 
S/N > 200. 

In the solar spectrum, we selected 54 unblended lines of Fe I and 18 lines of Fe II which cover 
a broad range of the line strength and excitation energy of the lower level. Despite the existence 
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of many sources of (7/— values for neutral iron, there is no single source that provides data for all 
the selected Fe I lines. We employ experimental (7/— values from Bard et al. (1991), Bard & Kock 
(1994), Blackweh et al. (1979; 1982a; 1982b), Fuhr et al. (1988), O'Brian et al. (1991). Multiple 
sources of 5/— value are available for each line of Fe II. We inspected 5 sets of data from Melendez & 
Barbuy (2009, hereafter MB09), Moity (1983, hereafter M83), Raassen & Uylings (1998, hereafter 
RU98), Schnabel et al. (2004, hereafter SSK04), and the VALD database (Kupka et al., 1999). Van 
der Waals broadening of the iron lines is accounted for using the most accurate data available from 
calculations of Anstee & O'Mara (1995), Barklem & O'Mara (1997), Barklem et al. (1998), and 
Barklem & Aspelund- Johansson (2005). 

The Sun. With the solar model atmosphere 5780/4.44, the non-LTE effects are small for the Fe I 
lines and negigible for the Fe II lines, such that the non-LTE abundance correction amounts, on 
average, Anlte = log snlte — log slte = +0.03 dex for Fe I and smaller than 0.01 dex for Fe II. A 
depth-independent microturbulence of 0.9kms~^ was adopted in the calculations. The iron abun- 
dance was derived from line profile fitting. For the Fe I lines, we obtained the mean non-LTE 
abundance logeFei = 7.56 it 0.09. The statistical error which is defined as the dispersion in the 
single line measurements about the mean, a = \/T,{x — x^j^ j (n — 1), is considerably too high. This 
is, most probably, due to the uncertaunty in (7/— values because the abundance scatter is largely 
removed in a line-by-line differential analysis of solar type stars as shown by Mashonkina et al. 
(2011). For example, the Fe I based abundance [Fe/H]/ = 0.11 ± 0.03 for /3 Vir (6060/4.11). We 
find that the average Fe II based abundance depends significantly on the used source of 5/— values: 
logepeii = 7.41 ± 0.11 (SSK04), 7.45 ± 0.07 (VALD), 7.47 ± 0.05 (MB09), and 7.56 ± 0.05 (RU98, 
M83). It is worth noting that the most recent laboratory g/— values of SSK04 lead to the highest 
abundance scatter. 

Large statistical errors and significant systematic discrepancies between different authors make 
the situation with oscillator strengths for the visible Fe I and Fe II lines to be unacceptable. Inde- 
pendent of either ID or 3D modelling, no firm conclusion can be drawn with respect to whether or 
not the Fe I/Fe II ionization equilibrium is fulfilled in the solar atmosphere and how realistic the 
non-LTE calculations for iron are. The astrophysical community shall ask atomic spectroscopists 
for new accurate measurements. 

With the cited 5/— values for Fe I and the data of RU98 for Fe II, the ionization equilibrium 
between Fe I and Fe II is matched consistently with the solar log g = 4.44. Is it possible to achieve 
a similar consistency for any other star using the same line list and the same atomic data? 

Procyon. Procyon is among the very few stars for which the whole set of fundamental stellar 
parameters except metallicity can be determined from (nearly) model-independent methods. Allende 
Prieto et al. (2002) obtained T^g = 6510±49 K and log g = 3.96±0.02. Microturbulence velocity 
of T4nic=1.6 kms~^ was determined in this study from the requirement that the non-LTE iron 
abundance derived from Fe I lines must not depend on the line strength. 

We find that the mean LTE abundance determined from the Fe I lines is 0.14 dex lower than 
that from the Fe II lines. The non-LTE abundances from the individual Fe I and Fe II lines are 
presented in Fig.[2l with the mean values of logepei = 7.55 ± 0.09 and logepeii = 7.62 it 0.06. Non- 
LTE removes, in part, a disparity between two ions and leads to consistent within the error bars 
abundances from two ionization stages. But the abundance error is as large as for the Sun, and this 
leaves a space for speculations on either one needs to improve further the non-LTE line formation 
modelling, or to revise upwards Tgfj (change in Tefj by 80 K fully removes a disparity between Fe I 
and Fe II), or to apply a new generation model atmosphere based on hydrodynamic calculations. 
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Figure 2: Non-LTE iron abundances from the Fe I (open circles) and Fe II (filled circles) lines in 
Procyon plotted as a function of Wots- The dashed line indicates the mean abundance derived from 
the Fe I lines and the shaded grey area its statistical error. For Fe II, gf— values are from Raassen 
and Uylings (1998). 



Table 1: Non-LTE abundance corrections (dex) for the selected lines of Fe I and Fe II. Multiplet 
numbers are indicated in the third string. 
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5. Departures from LTE depending on stellar parameters 

The non-LTE calculations were performed for the small grid of model atmospheres with stellar 
parameters characteristic of early F to late A-type stars: T^d = 6500 — 8500 K, log g = 3 and 4. 
In this stellar parameter range, the general behavior of the departure coefficients of the Fe I and 
Fe II levels is independent of effective temperature and surface gravity and very similar to that 
shown in Fig.[TJ Non-LTE leads to weakening the Fe I lines and to the opposite effect for Fe II. 
TablelU presents the non-LTE abundance corrections for the representative lines of Fe I and Fe II 
with various -Eexc- 

The non-LTE effects are only minor for Fe II, such that the abundance correction does not 
exceed 0.03 dex in absolute value. For Fe I in the log g = 4 models, ANLTE < 0.12 dex and, on 
average, decreases towards higher Tegf. For the lower gravity models, ANLTE < 0.20 dex and, on 
average, grows towards higher Tes- It is worth noting that the non-LTE abundance correction is 
substantially different for different lines, at any given temperature. 

Our theoretical results were compared with the non-LTE predictions of Rentzsch-Holm (1996) 
in common stellar parameter range, T^q = 7000 - 8500 K, log g = 4 and 3.5. There is no discrepancy 
for the Fe II lines, while the departures from LTE for Fe I are stronger in Rentzsch-Holm (1996) 
than in this study. For example, in the model 8500/4, Rentzsch-Holm (1996) gave the mean non- 
LTE correction of 0.17 dex which is 0.1 dex larger than our value. In contrast to our results, the 
non-LTE correction from Rentzsch-Holm (1996) grows towards higher Tefj independent of stellar 
surface gravity. In both studies, the departures from LTE grow towards lower gravity, however, 
with log g = 3, we obtain the smaller non-LTE corrections than Rentzsch-Holm (1996) did with 
log g = 3.5. The difference between two studies is in the use of different model atoms. Our model 
atom of Fe I is much more complete than that of Rentzsch-Holm (1996), where only 79 terms of 
Fe I are included and most high-excitation levels with -Eexc > 6 eV are absent. This explains why 
the overionization of Fe I is stronger in the Rentzsch-Holm (1996) calculations. 

6. Conclusions 

• Completeness of the model atom is important for correct calculation of the statistical equi- 
librium of iron. 

• In the stellar parameter range characteristic of middle F to middle A-type stars, LTE under- 
estimates the element abundance derived from Fe I lines, by 0.02 - 0.1 dex for dwarfs depending on 
the line and stellar temperature and by up to 0.20 dex for giants. 

• LTE is as good as non-LTE for Fe II. 

• The uncertainty in available 5/— values for visible lines of Fe I and Fe II is uncomfortably 
too high, such that solar and stellar iron abundances cannot be determined with the accuracy 
better than 0.1 dex. The astrophysical community shall ask atomic spectroscopists for new accurate 
measurements. 
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